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Background 


Since before man ventured into space, concern was expressed for the 
advent of an impact of the space vehicle by meteoretic particles. 
Then, the concern was for asteroidal particles or cometary debris. The 
average particle density was of the order of one-half gram per cubic 
centimeter (0,5 g/cc) and the average encounter velocity was on the 
order of twenty kilometers per second (20 km/sec). Early probes such 
as the Explorers and Pegasus led to laboratory experiments of impact 
into a variety of materials with a wide range of projectiles and 
velocities. 

In support of this early research, this experimenter proposed a pene- 
tration criteria for a variety of metal targets at the 7th Hyper- 
velocity Impact Symposium at Tampa in 1965. (1) This penetration 
criteria became known as the Fish-Summers penetration equation for 
thin targets. Further experimentation by others (2,3) substantiated 
the findings and led to adoption of the penetration criteria for 
single sheets into the Meteoroid Hazards Monograph (4). 

As our space effort grew, early distribution theories of meteoroid 
population were modified to reflect the expanding knowledge of our 
solar system. 

In recent years, however, the numbers of orbiting satellites, spent 
components, collisions and explosions has populated the near earth 
orbit with debris potentially more hazardous than the average meteo- 
roid debris. This new debris has an average density of aluminum (2.78 
g/cc) and and and overage encounter velocity of 10 km/sec. 

The space station will require many hours of work in this environment 
and there is concern over hazard to the assembly personnel. The 
proposed 'hard suit' design utilizes 6061-T6 Aluminum for most of its 
exposed area. The aluminum surface will be treated for thermal and 
radiation control. The basic thickness of this suit will be on the 
order of 1.78 mm (0.070 inches). 

Since the earlier work in the sixties was performed generally with 
2024-T4 and 1100-0 aluminum it was felt that some ballistics testing 
on 6061-T6 aluminum to verify its ballistic limit was warranted. This 
would allow updating the penetration criteria for the newer material. 
Further, it was felt that some impact testing of coated specimens was 
necessary to determine what if any debonding occurred during the 
impact process. 

Test Procedures: 

Preliminary Test: 

Before testing the newer suit materials, a shoulder section of current 
suit material was to be impacted. This operational section was tested 
in a configuration similar to that tested at Johnson Space Center at 
least as for as particle size was concerned. An aluminum sphere of 
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0.79 mm was selected to impact at approximately 5.5 km/sec. 

The section was sealed, pressure tested and monitored for average leak 
rate. Since the section could not be subjected to extremes in pres- 
sure, an arrangement was devised to allow monitoring pressure at all 
times from outside the evacuated range before shooting. This was 
accomplished by using an evacuated accumulator tank and double gauges. 
See figure (1). 

As the range was pumped down the closed suit section internal pressure 
began to rise. A valve on the external evacuated tank was opened to 
lower the suit pressure below its limit. When 3 to 5 mmhg range 
pressure was reached the accumulator tank valve was shut off to trap 
the suit pressure at 8 psig on the gauge located inside of the range. 
Under these conditions the round was fired. 

A small block of wood was placed inside of the suit section to catch 
any debris resulting from the impact. Rapid pressure loss followed 
the impact. The projectile was observed to have broken up into 3 or 4 
pieces which imbedded in the wood. The impact hole was readily seen 
but was easily patched with tape from the inside of the suit. 

( Comments on this initial test ) 

This suit shot was an effort to try and reproduce fabric shots in an 
untitled report from JSC of tests performed on various components of a 
candidate space suit. 

Several questions arose when reading this hazard assessment of a 
candidate 8psi space suit configuration. Some of these questions were 
answered in a phone conversation with Bert Cour-Palais. One of the 
questions posed was why the energy levels for penetration of 6061-T6 
Aluminum were calculated using the hypervelocity equation cited in the 
report, namely: 

(1) tp = K x d**(1.06) x (pp/pt)**(0.5) x (V/C)**(.667) x (BH)**(-.25) 


Where d = diameter, cm 

pp <= projectile density > 2.78 g/cc 

pt = Target density > 2.78 g/cc 

V = Impact Velocity, cm/sec > 10E05 cm/sec 

C = Sonic Velocity " > 5E05 " 

BH « Brinnel Hardness > 95 


K = material constant 1.75 x 5.24 -> 9.2 
tp = threshold perforation, in cm 

This equation was developed by Rockwell and also JPL some years ago. 
The penetration equation used K=5.24. JSC in this report uses a modi- 
fied K which is 1.75 times the K=5.24, or K=9.2. Some earlier reports 
on meteoroid bumpers uses a factor of 2.0 times for a K=10.48 for a 
bumper of zero penetration. 

Several other groups use a modified Fish-Summers equation, namely: 

(2) tp = K1 x (pp)**(.1667) x (mp)**(.352) x (Vp)**(.875) 
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Where pp = Projectile density g/cc > 2.78 g/cc 

mp = Projectile mass g > 

Vp = Impact Velocity km/sec > 10 km/sec 

K1 = Target Material Constant Alum > 0.54 to 0.57 

tp = threshold perforation in cm 

The 1969 Meteoroid Environment Model used K1 = 0.54 for Aluminum. A 
1984 MSFC report by Dallas Smith on the debris impact uses Kl=0.57. 

Plots of these equations are shown in figures (2) & (3). If we look at 
a plot of the F-S equation with Kl=.54 and overlay a plot of the JSC 
equation with K=9.2 this is shown in figure (3). These plots show 
Thickness of penetration versus projectile diameter for impacts at a 
constant of 10 km/sec. 

The foregoing discussion on which formulas are being used and the true 
ballistic limit of the 6061 T6 suit material will be determined in the 
plan to impact 6061 T6 Aluminum with Aluminum projectiles to determine 
the ballistic limit. This B.L. can then be used to determine the 
probability of no penetration of the proposed suit for the space 
station. The ballistic limit shots will preceed the shots to determine 
the effect of impact upon the protective coatings being considered.) 

Ballistic Limit of 6061-T6 Aluminums 

An initial series of impacts using various sizes of projectiles was 
conducted into target plates of 6061-T6 aluminum. The projectiles all 
of aluminum and spherical. Figures (4) & (5) show a cross-section of 
sheets impacted with 0.79 mm (0.03125inch) and 1.5^ mm (0.0625 inch) 
aluminum spheres at approximately 5.5 km/sec. However, the 

determination of velocity is very poor with projectiles in this size 
range. The detectors are extinction type, i.e., they look for an 
interruption of a bright field with the passage of the particle. When 
the particle size is small the beam interruption goes unnoticed being 
lost in-the-noise, and if a stop signal is given the range counters, 
one has no garrantee that the signal was due to a particle passage or 
not . 

Therefore, it was decided to use 3.175 mm (0.125 inch) diameter 
aluminum projectiles since reliable velocity determination could be 
made. If we assume that the behavior of 6061-T6 aluminum to be similar 
to that of other aluminums we can use the Fish-Summers equation to 
predict the approximate perforation velocity limits and reduce the 
, number of shots required for this series of tests. 

If we take the modified Fish-Summers equation cited above, namely: 

(2) tp = K1 x (pp)«*(.1667) x (mp)**(. 352 ) x (Vp)**(.875) 


Where pp = Projectile density g/cc > 2.78 g/cc 

mp = Projectile mass g > 

Vp = Impact Velocity km/sec > 10 km/sec 

K1 = Target Material Constant Alum > 0.54 to 0.57 

tp = Thickness penetrated in cm 
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We may use either K1 = 0.54 or K1 = 0.57 initially to get an impact 
velocity near that assumed for penetration of the specimens chosen for 
the impact tests. Figure (6) is a plot of penetration thickness 
versus impact velocity for a one-eighth inch diameter aluminum pro- 
jectile into an aluminum target with K1 = 0.54. We see that a 
projectile of this size with a velocity in the range of 5 to 6 km/sec 
should penetrate on aluminum target ranging from 0.9 cm to 1.0 cm. 
Using 0.953 cm (0.375 inch) 6061-T6 aluminum plate as a target 
material to match the 0.3175 cm (0.125 inch) aluminum projectile 
should allow determining threshold perforation within the capabilities 
of the Vertical Gun Range Facility. 

Several rounds were required to zero in on the threshold Derforation 
limit, which was found to be 5.02 km/sec. This suggested that the 
value for K1 target constant lor 6061-T6 aluminum should be 0.57. 
Figure (7) is a plot of thickness versus impact velocity for 1/8 inch 
diameter aluminum projectile and a target constant K1 = 0.57, The 
threshold perforation point is noted on the plot. 

Figure (8) is a plot of suit thickness versus debris diameter with an 
impact velocity of 10 km/sec. This suggests that for normal impact the 
suit thickness of 1.78 mm (0.070 inch) would be perforated iy d.S'bris 
diameters greater than 0.35 mm at 10 km/sec. 

If we consider the appearance of the target cross-section at 
threshold perforation, earlier work (a) pointed to the difference be- 
tween highly ductile and brittle materials. Consider figure (9), (a) & 
(b), where (a) is a drawing of a brittle material and (b) is that of a 
ductile material of the same density. 

We see that the brittle material, while still retaining a roughly 
hemispherical crater at threshold perforation, fails through cracks to 
the rear surface where spall material has been ejected by the 
reflected shock wave interaction. 

For the highly ductile material, (9b), the crater shape has distended 
to a conical section, failure occuring through a hole at the base of 
the crater, and the spall products retained in petals below the 
crater. 

If we compare these figures with figure (10) of the threshold 
perforation of 6061-T6 aluminum, we see some distension and spall 
retension like a ductile material but with some cracking through like 
that of a brittle material. This is indicative of behavior of 6061-T6 
aluminum with a balance between brittleness and ductility. 

What happens at other than normal impact? Certainly all impacts will 
not be normal to the surface. Figure (11) shows 6061-T6 aluminum 
impacted at 30 degrees to normal. The point is just below threshold 
perforation. Compare this figure with that of the previous shown fig- 
ure (10) for threshold perforation at normal impact. The cross- 
sections are almost identical in appearance. The velocity is about 0.5 
km/sec higher than the earlier impact. The cratering phenomena occurs 
due to the pressures developed at impact which sets off shock waves at 
right angles to the surface at the impact point. For all but extremely 
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low angles of Impact, the craters will not betray the impact angle. 

The pressures developed will be governed by the normal component of 
the particle velocity, that is, a function of the cosine of the impact 
angle. We could probably rewrite the modified Fish-Summers to include 
this angular dependence. 

(3) tp = K1 x (pp)**( .1667) x (mp)**(.352) x (V)**(.875) x (cos(phi)) 

Where tp » Thickness penetrated 

pp * Particle density AL = 2.78 g/cc 

mp = ” mass 

V= " velocity 

phi = angle of impact (from normal) 

K1 o Target constant (0.57 for 6061-T6 AL) 

Impact of Coated Specimens: 

Determination of the best available thermal and radiation control 
coatings was arrived at by a separate study. Two candidate coatings in 
different configurations were presented for impact study. These were 
Aluminum and Gold in both polished and satin finishes. Also satin 
finished specimens of each were subjected to low-speed> high-mass 
impacts from tumbling with other materials such as tools, etc. These 
latter specimens had a large number of burnished areas resulting from 
these impacts. The specimens were dimensioned about 2 inch by ^ inch 
by 1/U inch thick. It was felt that maximum effect would occur if the 
high speed impacts were with particles of insufficient mass and velocity 
to penetrate the specimens but rather provide a reasonable crater and 
be as uniform as possible in the impact energies. 

With the 1 M inch thick plates of 6061-T6 aluminum coated specimens it 
was decided to use 1/16 inch (1.59 mm) aluminum projectiles at around 
5 km/sec velocity. This would keep the impacts below penetration point 
and not give rise to any large spall areas on their rear surfaces. 

Figures (12) to (17) show the results of these shots. Figure (12) is 
an overall view of all the specimens, the balance of the figures show 
the craters close up. The areas of highest stress in the coatings, 
right adjacent to the craters showed many stress lines in the coatings 
but no evidence of any debonding or cracking of the coatings. Whatever 
coatings are used it appears that all tested very well under high 
speed impact. 

Conclusions: 

The selection of 6061-T6 Aluminum for space suits for use on the space 
station would appear to be worthwhile. The relatively ductile behavior 
of 6061-T6 aluminum is better than a choice of a more brittle 
material. 

The coatings selected show no signs of failure occuring around the 
impact craters as other materials might v/ell do. 

The suit for long term orbit use will likely include insulation mater- 
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ials. These can only act beneficially in reducing ballistic limits. 
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Figure 17 

Polished Finish Gold Coating 
1/16" A1 Sphere Impact 
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10 short program to calculate penetratiun thicl.ness of 
?W aluminum space suits using threshold perforation formula 
30 author diet fish, nasa ames research center 19134 
40 Get the data if different from del aul Ls: (FOR ALUM) 

50 L)I = .0625 : HO = 2.78 : I I - .5/ 

60 DIM X(101> , Y ( 101 ) 

62 COMMON X ( ) , Y ( > 

65 CLS 

70 PRINT “This program calculates the Threshold per for at i on" 

80 PRINT "for debris impact into structures in near earth orbit" 

90 PRINT "Input will be asked for , a N' answer will select" 

100 PRINT "default values...." : PRINI : PRINT 

110 INPUT "Change diameter of projectile? (dlt 1/16 inch) CY/143 " 5 

120 IF A$ ="Y“ OR A$ — " y " THEN 140 

130 IF A$ = "N" OR At ="n" THEN 150 ELLE 110 

140 INPUT "Enter new diameter in decimal inches ( .XXX ) ";DI 
150 INPUT "Change projectile density (dlt ALUM— 2. 70) CY/NI " ; Bt 
160 IF B* =“Y" OR B* ="y" THEN 180 

170 IF Bt ="N" OR B* = "n" THEN 190 ELSE 150 

130 INPUT "Enter new projectile density in qrns/cc ( X.XX ) " ; K0 

1B5 Kl*=” ALUMINUM" 

190 INPUT "Change target constant Kl? (dlt is At UMINUM— .57) IY/H1 ; CF 

200 IF C$ ="Y" OR Ct «="y" THEN 220 

210 IF C$ ="N“ OR Ct — "n" THEN 230 ELSE 190 

220 INPUT "Enter new target constant temper ical data) (.XX) " ;l',l 

225 LINE INPUT "What is the new material? ";K1* 

230 D = D I * 2 . 54 : M — (l) 3) * ( . 5236) *RO : Ml — M (. 352) : 1)1 — RU (.16o7) 

2.40 DEF FNA ( X ) = K1 «M1 *D1 * ( X ( . 0/5) ) 

245 print 

250 PRIM! "tilling ar i ay with thickness and velocity values" 

260 FOR I - 0 to 100 
270 X ( 1 ) = 1/10 

2 00 Yd) = FNA ( X ( I ) ) 

270 NEXT I 

300 CHAIN " PLOT XY" ,5000 
997 END 


oRiGnys. 'is) 

OF POOR QUALFW 



C>M^JX V | 


10 'short program to calculate penetration thickness of 
10 aluminum space suits using threshold perforation formula 
30 'author dick fish, nasa ames research cenLer 1985 
40 ‘Get the data if different from defaults: (FUR ALUM) 

50 V = 10 : RD = 2.78 : K1 = .57 :K1 t=" ALUMINUM" 

60 DIM X (101 ) , Y ( 101 ) 

62 COMMON X ( ) , Y ( ) 

65 PR I NT sFRINT :FR1NT 

70 PRINT "This program calculates the Threshold perforation" 

O0 PRINT "for - debris impact into structures in near earth orbit" 

90 PRINT "Input will be asked for, a 'll answer will select" 

100 PRINT "default values " : PRINT : PRINT 

110 INPUT "Change velocity of projectile? (dft 10 km/sec) LY/NJ ";AT 

120 IF A* ="Y" OR AT =”y" THEN 140 

130 IF AT = "N" OR AT = "n" THEN 150 ELSE 110 

140 INPU1 "Enter - new velocity in decimal km/sec ( X.X ) " ; V 

150 INPUT "Change projectile density (dft 2.70) 1IY/N3 " ; BT 
160 IF ="Y" OR BT ="y" THEN 1S0 

170 IF BT = "N" OR BT ="n" 1 1 IEN 190 ELSE 150 

180 1NPU1 "Enter new projectile density in qms/cc ( X.XX ) " ; F<0 

190 INPUT "Change target constant Kl? (dft is ALUM 1 NUM= . 57 ) CY/N3 " ; CT 

200 IF C$ =" Y" OR CT =“y" THEN 220 

210 IF C$ = "N" OR CT ="n" THEN 230 ELSE 190 

220 INPUT "Enter new target constant (emperical data) (.XX) " ; Kl 

225 LINE INPUT "What is the new material? " ;K1$ 

230 VI = V' (.875) s M = (.5236)*R0 : D1 = ROM. 1667) 

240 DEF FNA ( X > = Kl *D1 * VI * ( X''. 352) 

245 PRINTS PRINT 

250 PRINT "Fi 1 1 ing array with target thickness vs projectile diameter values 

260 FDR 1 = 0 to 100 

270 X(l> = 1/200 

275 Ml = <X<1)'3)*I1 

280 Y ( I > = FNA (Ml > 

290 NEXT I 

300 CHAIN "PLOTXY" ,5000 
999 END 





ORIGINAL PAL':" 13 
OF POOR QUAL5Y Y 


10 short progi am Id calculate pel iet r a t. i on thickness of 
i\A aluminum space suits using Ber t Cuur -t'alais current f oriiiula 
30 author diet fish, nasa ames research center 1985 
40 Get the data if different from clef aul ts: (FOR ALUM) 

50 V = 10 s RO = 2. 78 : K1 - 5.24 : K1 1= “ ALUi II NUN" 

55 C = 5.1 : DU = 95. : PT = 2.70 
60 DIM X<101> , Y ( 101 > 

62 COMMON X ( > , Y ( ) 

65 PRINT :PF<INT SPRINT 

70 PRINT "This program calculates the Threshold per f or at i on" 

80 PRINT “for debris impact into structures in near earth orbit" 

90 PRINT " Input will be asked for, a ’N* answer will select" 

100 PRINT "default values...." : PRINT : PRINT 

110 INPUT "Change velocity of projectile- (dft 10 kin/sec) C i /N1 " iAT 

120 IF AT — " Y " OR AT ="y " THEN 140 

130 IF AT = "N" OR AT = "n" THEN 150 ELSE 110 

140 INPUT “Enter new velocity in decimal km/ sec ( X.X ) " ;V 

150 INPUT "Change projectile density (dft 2.7B) CY/N3 " ; BT 
160 IF BT ="Y" OR BT = "y" THEN 100 

170 IF BT = "N" OR BT = "n" THEN 190 ELSE 150 

100 INPUT "Enter new projectile density in gms/cc ( X.XX ) “ ; f<0 

190 I NPUT "Change target constant Kl? (dft is ALUM1NUM=5. 24) CY/N1 " ;CT 

200 IF CT ="Y" OR CT - " y " THEN 220 

210 IF CT = "N" OR CT = "n" THEN 230 ELSE 190 

220 INPUT "Enter new target constant temper i cal data) (.XX) ";K1 

225 LINE INPUT "What is Hie new material? ";K1T 

230 VI = (V/C) (.667) s B1 = (BH) -(.25) s D1 = (RO/PT) ' ( . 5) 

240 DEF FNA(X) - K 1 » L>1 * VI *B1 * ( X 1 . 06) 

245 PRINTS PRINt 

250 PRIN1 "Fi 1 1 i ng array with target tin cl. ness vs projectile cli aireier values" 
260 FUR 1 = 0 to 100 
2/0 X(I> = I/2U0 
275 Ml = X(l> 

280 Y ( I ) = FNA(Ml) 

290 NEXT I 

300 CHAIN “PLOTXY** ,5000 
999 END 
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REM* SLOT PROGRAM FOR 

REM* BY 

REM* Dick Fi Eh - 9/26/84 

REM* 

REM******-* ****** ****************************************** ******* 
REM 

REM USE THE FOLLOWING SPACE TO CALCULATE A 101 BY 101 

REM DATA ARRAY IN X(I) AND Y(I>. 

REM ENTER PROGRAM E<ELDW THIS SECTION IF ARRAY EXISTS. 

REM 

REM SAMPLE ARRAY— > CONCENTRIC El-LIPSES 

REM 

DIM Y (101 > ,X (101 ) :REM MOV DIM 10 CALLING PROGRAM 

REM 

FOR 1=0 TO 50 

Y ( I > =SIN( 1*6. 202/50) : X ( 1 ) =C0‘3 ( 1*6. 282/50) 

NEXT I 

REM 

FOR 1=51 TO 100 
J=I~50 

Y ( I) =. 5 xSINU *6.282/50) : X ( I ) = . 5*C0S ( J * 6. 282/50) 

NEXT 1 

REM 

REM *** CALLING PROGRAM MUST DIM XO & Y() & COMMON *** 

REM***** ********************************** ** *********** ********** 
REM* ENTER TITLE OF PLOT AND AXIS LEGENDS * 

R[rjvj* *X* X**XXttX**XX XX ****** **XXi(XXK*X*ftX**XXX**X ft*** ***Xftft*X ****** 

REM WV ENTER PLOT HERE AS SIJB 10 MAIN MATH PROGRAM WV 
REM 

XT I TLEP=SPACEP (50) : YTITLEP=GPACEP (20) : FT I TLEP=SPACEP (50) 

LINE INPUT "ENTER TITLE OF PLOT (50 CHARACTERS OR LESS) " ; PP 
LINE INPUT "ENTER Y-AXIS TITLE (20 CHARACTERS OR LESS) " ; YP 
LINE INPUT "ENTER X AXIS TITLE (50 CHARACTERS 
MIDP (Y I I VEEP, (LLN(YTIILEP) --LEN(YP) ) /2>=YP:REM 
MIDP ( X 1 1 TLEP , (LEN (XTITLEP) -LEN (XP) ) /2)=XP:REM 
M I DP < P T I TLEP- , (LEN (PT ITLEP) - LEN ( PP ) ) / 2) —F $ : REM 
INPUT "DO YOU WANT A(UTO) OR M(ANOAL) SCALING 


OR LESS) " ;XP 
i THIS SECT I UN 
: SCALES H IE 

: TITLE SPACE 

" } NP 


IF 

IF 


NP= "M" 
NT -"A" 


OR 

OR 


in 


NP = 
Ni=“a" 


THEN 5110 
THEN GOLUB 


GOTO ‘ 
I NPLJ r 
INPU 1 
INPUT 
INPUT 


1 2 05 
“ENTER 
"{ NTLR 
"ENTER 
"ENTER 


MAXIMUM VALUE OF 
MINIMUM VALUE OP 
MAXIMUM VALUE OF 
MINIMUM VALUE OF 


6330 ELSE 5070 

Y* AXIS"; YMAX 
Y-AXIS"; YNIN 
X AXIS" ; XMAX 
X- AXIS";XMIN 


REM 

REM* ***** ****#« * #*****#« ****** ************************** ********* 
REM* MAIN BODY OF PROGRAM SI ARTS HERE * 

F<EM* ****** ******************************** ******* ** ******** x *.+■ *** 


5190 REM 
5200 REM 

5205 CLS: LOCATE 1,19 

5210 PRINT PTITLEPsREM PRINT TITLE OF PL01 

5260 FOR LNNO=0 TO 180 STEP 18 

5270 YN=YHINt (( (YMAX- YMIN) / 180) *( 1B0LNNO) ) : REM NORMALIZE Y 

5200 YD=ADD( (YMAX- YNIN) / 180) : REH VALUE OF EACH LIME 

J290 GOSUB 5930: RLM PRINT VERTICAL SCALE 

5300 NEXT LNNO 

5305 GDSUB 6040: F(EM PRINT Y TITLE 

5310 GOSUB 50 1U: REIT PLOT DATA 



! 1 ill 


11) 1 ijl'H'j : i n:x I 


l iil.j'i 
; Hi X I J 


JCi .1.1- 


5370 130SUB 6130: REM 

5450 FOR J = l TO 20: I UR K=1 

545b SVG TEH 

5460 REM - — 

54 70 REM* *** ** ******************** K ***** ************ *'**>.**fc*>»* l 'fc*''- : **** 
J4S0 REM* END OF MAIN BODY- SUBROU I IMEG START HERE * 

5490 REM** * *************************** **************** *** »* ********** * 


ORIGf^L pf.VZ 1‘3 
OF POOR QUALJVY 


5500 REM 

5770 REM* ********************************* ** K * l< '** ********************* 
5780 REM* PLOT DAI A GUDROU TINE * 

5790 F\EM** ************************** K ****** x a ***** *' 1 ‘‘ ********* * 14 * * ***** 
5800 REM 

Y5= 1 3+ 1 80 - FIX ( (V (0) - YMII'I) * 180/ (YMAX-YMIN) ) 

XS=60*F I X ( (X(0) XMIN) *560/ (XNAX-XMIN) ) 

PSET (XB,YS),7 
FOR 1=1 TO 100 

YP=j 1 3+ 1 80-F 1 X C (Y(I) -YMIN) * IBB/ (YMAX-YMIN) ) 

XP=60+F I X ( ( X ( I ) — XMIN) *560/ ( XMAX— XMIN) ) 

IF YPC13 OR XP>620 THEN 5870 
LINE — ( XP , YP ) 

NEXT I 

LINE (60, 13) - (620, 193) , ,B 
RETURN 

5080 REM 

5090 REM*** ******* * ************************************* ****** ******** 
=,900 REM* PLOT Y SCALE * 

59 1 0 REM*** ******************* ******* ******* ************************* R 
5920 REM 


5930 

RT=" R20 " : LT= " L 20 “ 

5940 

LOCATE LNNQ/ 10*2+2,2 


PRINT USING "THMT. TNI ” 5 YN 

59 612 

1 >SET ( 60 , LNNOH 13 ) , 7 : DFtAW F< T 

59/0 

PSET (620,LNNO* 13) ,7s DRAW L4> 

5980 

RET URN 


5810 

5820 

5830 

5832 

5834 

5836 

5837 
5850 
5B60 
5870 
5875 


5990 REM 

<_.O00 REM* ***************** 
6010 REM* 

6020 REM*** ******* « ** ** *** 


*********** k'- ******* K ***************** R * * * # 
PRINT Y-AXIS TITLE * 

**************** ***** ** * *************** * * *** 


6030 REM 


6040 LOCATE 3 , 1 

6050 FOR 1=1 TO LEN(YTITLET) 

6060 CHART -MI DT (YT 11 LET , I , 1 > 


6065 PRINT CHART 

6070 NEXT I 

6080 RETURN 

6090 REM 

6100 REM***************************** ** ********* *********** ** 
6110 REM* PRINT X SCALE SUBROUTINE 

6120 REM*** ******** ** * * k * ***** ******************** * * ** * ** * * * +: 


* *** ** *** 
* 

« k k * « >: x * * 


6130 REM 

6140 UT = " U 1 0 " : DT= " D 1 0 " 

6150 LOCATE 23,4 

6160 FOR 1=0 TO 10 

6170 PRINT USING ###" j ( XMAX- XMIN) * I / 10 «• XMIN; 

6172 NEXT I 

6174 FOR I-l TO 10 

j 1 7 5 l ‘SE J ( 56 * I i 4 , 1 '7 3 ) , / s DT \ ATJ UT 

6 1 76 PSE I (56* 1*4,13) ,7 : DRAW DT 

6180 NEXT I 

6200 REM 



62 1 0 RE M » * If i. ;; * r: x i: k a i: a „ r. f •». K >. >' * * * !•: >. a K • > > ■ > « - * ; ’• '• **** * * * 5 ;; • Y " - * •• 

62 2D REM* PRINT X AXIS Til LI. * 

6230 REM*** * a* a a a k«*m*k a **************** **■*****■ ****** *********** *'■»*** 


62A0 REM 

6250 LOCATE 24,19 

^260 PRINT XT ITLE$; 

6265 LOCATE 1 , 1 

6270 RETURN 

6280 REM 

6290 REM***** ***************** *******“ ********** ********* *********** * ** 
6300 REM* AUTO SCALING SUBROUTINES * 

6310 REM*** *** **************** ********************************** ******* 
6320 REM 

6330 YHAX=Y (0) s XMAX = X (0) : YMIN-Y (0> : XM1N-X <Q> 

6340 FOR 1=1 TO 100: REM FIND XMAX & YMAX 

6350 IF Y ( I > >YMAX TMEU YMAX=Y(I) 

6360 IF YlIXYMIN THEN YMIN=Y(I> 

6370 IF X ( I ) >XMAX TMEN XMAX=X(I> 

6380 IF XdXXMIN THEN XMIN=X(I> 

6390 NEXT 1 

6400 RESTORE 6610: REM 

6410 HSD= (YHAX-YMIN) / 10:REM SCALE THE Y AXIS 


6420 

6430 

6440 

6450 

6460 

6470 

6480 

6490 

>500 

6510 


FOR L=-2 TO 4 

FOR K*1 TO 3: READ J 

IF MSD<=J*10"- (1> THEN MSD= J GOT O 6460 

NEXT Ks RESTORE 6610: NEXT I 
FOR 1=10 TO -10 STEP -1 

IF (YMAXOHMSD) * (YMAX>I*HSD-. 79999*MSD> TMEN YMAX=I*MSD 
NEXT I 

YM I N= YMAX -- 10“ M'SD 

RESTORE 66 10: REM- 

MSD= (XMAX— XMIN) / 10: REM SCALE THE X AXIS 


6520 

6530 

6540 

6550 

6560 

6570 

6530 

6570 

6600 

6610 


FOR 1=2 TO 4 

FOR K=1 TO 3: READ J 

IF HBD<=J*10'I TMEN M f SD=J * 1 U I : GOTO 6560 
NEXT K: RESTORE 6610: NEXT I 
FOR 1 = 10 TO - 10 STEP -1 

IF (XMAX< = I*MSD) * (XMAX>I*MSR- . 9V999*M3D) TMEN XMAX-I+MSD 
NEXT I 

XMIN=XMAX-10*N3D 
RET UFtN 

DATA 1 ,2,5: REM CHANGE THIS DATA LINE If 


6750 REM 

6760 REM* •K-***##K*****#*tth#*HKKa«K«**«*K#*K-tCKI' i -K#K****-K** +: #K***K* , !*' ; ' 11 * K * 
6770 REM* LIST OF VARIABLE NAMES USED IN THE PROGRAM * 

6780 REM*** **** a************ ********* * * ********* ***** ******* *** * * K * * 

Y AXIS DATA AF<RAY 
X AXIS DATA ARRAY 
MAXIMUM VALUE OF Y AXIS 
MINIMUM VALUE OF Y AXIS 
MAXIMUM VALUE OF X AXIS 
MINIMUM VALUE OF X AXIS 


6790 

6800 

REM 

REM 

Y ( I > 

6S10 

REIT 

X(I) 

6820 

REM 

YMAX 

6830 

REM 

YMIN 

6840 

REM 

XMAX 

6050 

REM 

XMIN 


6860 REM 
6870 REM 
68B0 REM 
j890 REM 
6900 REM 
6910 REM 
6920 REM 


Lt "NO 

RSTOf 

LST ACE 

YN 

YD 

MSD 

J 


Y AXIS LINE NUMBER 
CHARACTER TAD STOP 
PRINTER LINEFEED IN U/72 INCH 
NORMALIZED Y AXIS VALUE 
DELTA VALUE EACH Y A >. I G 1 INC 
MINIMUM SCALE DELTA 
SCALE DELTA 



6/30 

IT il 

r; 

69 TO 

REIT 

i 

6'/. 60 

REIT 

IT 

6 >L0 

REM 

FT 1 1 L E 1 

69 /B 

REM 

Y* 

j,700 

REM 

YI11 LET 

6990 

REM 

x* 

7000 

REM 

X IITLL* 

7010 

REM 

CHAR* 

7030 

REM 

N* 

7060 

9510 

REM 

REM* 

I , .1 , K , L , ?<M 


INDEX lT'IIRER 
INDEX COUNT nr< 

FLO! TITLE 

CENTERED PLOT TITLE STRING 
V AXIS TITLE 

CENTERED Y-AXIS TITLE STRING 
X AXIS TITLE 

CENTERED X-AXIS TITLE SIRING 
CHARACTER PRINTED AT RSI OF 
INPUT VARIABLE FOR SCALING 
! CLOT OF DATA 


ORIGINAL 

OF POOR QUALITY 


